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Abstract 

The radiative leptonic Be — > decay is analysed in context of 2HDM. It is 

shown that with large values of tan(5, the contributions of Model 11 to the decay rate 
exceeds considerably the Standard Model ones, while the contributions of Model 1 
overlap with the Standard Model predictions. 
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1 Introduction 



The observation of the i?c-nieson by the CDF collaboration in the channel Be J/ ipiv, 
with ground state mass = 6.4 ± 0.39 ± 0.13 GeV, and lifetime r{Be) = 0.46lg;i^±0.03 
ps [1], has stimulated up the investigation of the properties of Sc-mesons theoretically, and 
experimentally on a new footing. The particular interest of this observation is related to 
the fact that, the meson ground-state with hc{hc) can decay only weakly, thus providing 
the rather unique opportunity of investigating weak decays in a heavy quarkonium system. 
Moreover, studying this meson could offer an unique probe to check the perturbative QCD 
predictions more precisely, and one can get essential new information about the confinement 
scale inside hadrons. The theoretical study of the pure-leptonic decays of i^c-meson, such as 
Be — > Ivf, {i — e, II, t) can be used to determine the leptonic decay constant /b^ [2], as well as 
the fundamental parameters in the Standard Model (SM), such as the Cabibbo- Kobayashi- 
Maskawa (CKM) matrix elements which are poorly known at present. Nevertheless, the 
well-known " helicity suppression" effect make an experimental difficulty in the measurement 
of purely leptonic decays of Be- Although, the Be — > ti7^ channel is free of the helicity 
suppression, the observation of this decay is experimentally difficult due to the efficiency 
problem for detecting the r lepton. 

Recently, the radiative leptonic Be — > ii^ej decays {i — e,ii,T), received considerable 
attention as a testing ground of SM and "new physics" , where no helicity suppression exists 
any more [3-5]. 

Among various radiative leptonic decays, the Be — >■ tz/t-7 decay provokes special inter- 
est, since the SM predictions has been exploited to establish a bound on the branching 
ratio of the above mentioned decay of order 10~^ [5,6], and therefore can be potentially 
measurable in the up coming LHC B-factories, where the number of Sc-mesons that will 
be produced are estimated to be fa 2.0 x 10^^ [7,8]. This will provide an alternative way to 
determine the decay constant fs^ and the CKM matrix elements. The decay Be — > ti7^7 
receive two types of contributions: internal bremsstrahlung (IB), and structure-dependent 
(SD) parts. The IB contributions are still helicity suppressed, while the SD ones contain the 
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electromagnetic coupling constant a but they are free of helicity suppression. Therefore, 
the radiative decay rates of Be — > (.v^^y could have an enhancement with respect to the 
purely leptonic modes of Be — > Ive, due to the SD contributions, thus it enable to establish 
"new physics" beyond the standard model. 

In this work, we will study the radiative leptonic B,, — > t^v7 decay in the framework 
of the two-Higgs doublet model (2HDM) [9-11] at large tan(3. The so-called Model I and 
Model II are considered, which are differ only in the couplings of the charged Higgs bosons 
to fermions. Subsequently, this paper is organized as follows: in section 2, the theoretical 
formalism relevance for the Be tz7^7 decay in 2HDM is presented. Section 3, is devoted 
to the numerical analysis and the discussion of the results. 

2 Formalism for the tiJ^-j decay 

In the Standard Model (SM), the decay Be — > t^v7 can be studied to a very good 
approximation in terms of four-fermion interactions. The effective Hamiltonian relevant to 
the process Be —>■ ti7t is: 

G 

Heff = -^ai1/c6C7^(l-75)&T7^(l-75)'^r, (1) 

where Gp is the Fermi constant, Veb is the CKM mixing element, ai is a QCD corrected 
factor, which is equal ai ~ 1.13. However, in the next discussion we will put ai ~ 1 

The emission of a real photon in leptonic decays of heavy mesons Be ti^tI can 
proceed via the two mechanisms mentioned in Sec. I. For the IB amplitude, the charged 
Sc-meson emits leptons via the axial- vector current, and the photon is radiated from the 
external charged particles. On the other hand the SD amplitude is governed by the vector 
and axial vector form factors, in which the photon is emitted from intermediate states. 
Gauge invariance leaves only two form factors /i,2(p^) undetermined in the SD part. The 
possible diagrams for this two mechanisms are shown in Fig. 1. Following this framework, 
the general form of the gauge invariant amplitude corresponding to Fig. 1 can be written 
as [6]: 

M{Be^Tl7r-f)^Mi + M2, (2) 
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where Mi and M2 represent the contributions of "inner bremsstrahlung" (IB), and "structure- 
dependent" (SD) parts, given by: 

Ml = le-^ VcbfBcmrSa u{pi) < W 1 - 75) v{p2) , (3) 



where Sa is the photon polarization vector, pi, p2, and q are the four momenta of r, i/^, 
and 7, respectively, /s^ is the Sc- meson leptonic decay constant, /i,2(p^) corresponding to 
parity conserving and a parity violating formfactors, p = Pb^ = Pi + P2 is the momentum 
transfer to lepton pair. The necessary matrix elements related to the Jb^, and to the hadron 
transition form factors /i,2(p^) are defined as follows [4]: 

(0 I 07^756 I B,) = -i/b^Pb^. , (5) 

(7(?) I c-fab I Bc{p + q)) = e^^^^eap^e^pPq^ , (6) 

fJp^) 

iiiq) I C7a75& I -Bc(p + g)) = -«e — —Sa [ga/sip-q) - Paq/s] ■ (7) 

We want now to consider the ru^'-f decay in the context of a 2HDM with no flavor 

changing neutral currents (FCNC) allowed at the tree level, i.e. Model I and Model II. The 
interaction lagrangian of fermions with the charged Higgs fields in both models is given by 
[12]: 

L = 2v?MvK I^VijiJiuiXuiil - ^5)dj + VijmdYui{l + 75)^^- 

+ meZP{l + -f5)?jH^ + h.c., (8) 

where gw is the weak coupling constant, is the W- boson mass, H"^ is the charged 
physical field, and Vij is the relevant elements of CKM matrix. In model I, X — cotp, 
Y — Z — —cotp and in model II, X — cotp, and Y — Z — tanp. 



The decay tVtI in 2HDM proceeds through the same Feynman diagrams (which 

are displayed in Fig. 1) that mediate the process in SM, except the VF-boson is replaced 
by the charged scalar Higgs boson H^, i.e. {W H^)- 



Fig. 1 

The matrix element corresponding to the diagram {W — > if^) where the photon is 
radiated from r lepton, (see fig.l) is: 

»fe){^?^i|±f^}(l-T5).fe), (9) 

where is the leptonic decay constant of Be meson, defined as: 
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(0 I CTsfe I B,{p + q)) = -ih^ (10) 

{nib + rric) 

While the contribution of the structure dependent part to the B^. tPtI decay, i.e., 
when photon is radiated from initial quark lines, due to the charged Higgs exchange can be 



obtained by considering the following correlation function: 

V2' 



< I T 



c(0) (m,Y{l + 75) + m,X{l - 75) )K0) 



rj{x) I i?e > X 



«(Pl)7/3(l -75)^^(P2), (11) 

where J^i^) is the electromagnetic current for 6 or c quarks. The hadronic matrix elements 
involving the scalar and pseudoscalar currents in Eq.(ll) are parameterized such that: 

/ dW^^ < I T{-c{0)^,b{0)rJ{x)) I B^ip + q) >= fn^ ""f- (12) 
J [mb + mc)p.q 

I d^xe'^^ < I r(c(0) b{0)J^\x)) I Se(p + g) >= 0. (13) 

The parameterization of the hadronic matrix elements given in Eqs. (12,13) are particularly 
well suited for our purposes since in the 2HDM, the vertex 6(c) ~ (1 —75) or 6(c) ~ (I + 75), 
hence, the vector part of this correlator is zero, and the active part of this correlator is the 
axial-part given by: 

^ -^'71 ^"^^^ "^i^l (m, + mj ^ 

^^(Pi)^(l-75)^(P2), (14) 

and the total matrix element in 2HDM becomes: 

V2 Mfj {nib + nic) 

At this accuracy it is easy to check that the modified total amplitude for the radiative 
leptonic B-decays of B^. — > tUtJ is gauge invariant: 

M(^totai) {B, ^ TJ7,7) = + M2, (16) 

where 

1 = VcbfEclTlrSa C u{pi)\ ^1 " 75) v{p2) ■ 17) 

V2 I 2pi.g p-q] 



Therefore, in this model the charged Higgs contribution modifies only the so-called Mi part 
of the SM, and it does not induce any new contribution to the so-called M2 (see Eq.3): 

The 2HDM is sensitive to two basic free parameters, namely tan(3, and the charged 
Higgs mass Mh- If we formally set Z — > in Eq. (18), the resulting expression is expected 
to coincide with the Be rurj decay, which was investigated in the framework of SM [6]. 

After lengthy, but straightforward calculation for the squared matrix element, we get: 

I Mtotai I' = I Mr I' + 2 Re [Mr^M^l + \ M2 \\ (19) 



where 



I' = — I V;6|V(-4/|m2)|C2^^^'' 



2 " {pi(iY{p<iY 

X ^p'ipiP2){piqf + ipqf [(piP2)(2m^ - p^q) + Mim^ - 2piq) 
+ {pq)ipiq) [ipP2)ipiq) - (pPi)(4piP2 +P2?)] \ , (20) 



2 Re 



X 



— I 
2 ' 

if2ip' 



m 



V;.|V(-16/^,m^)C2^^^ 
P^{Piq){p2q) + {pqf 



Be. 



{piq){pq) 
f2{p') 



2 {PiP2 + P2q) - '^^f-^ {P2q) 

IT^Bc ^Bc 



[{PP2){Piq) + (PP1)(P2?)] 



(21) 



Mo 



Vcb I'e^ie 



flip 



ml 



2^ |2 



+ 



f2{p 



ml 



2^ |2- 



X {(pP2)(pq)(piq) + (P2q) [(ppi)(pq) -p^(pig)]} 



(22) 



All calculations have been performed in the rest frame of the Be meson. The dot products of 
the four-vectors are defined if the photon and neutrino (or electron) energies are specified. 



The Dalitz boundary for the photon energy Ej and neutrino energy E2 are defined as 
follows: 



i2 ^2 



mn — mz. , , 

< R < — ^2 ^ . 23 



It is now straightforward to work out the expression for the differential decay rate in the 
lepton and photon energies: 

= 64^^' I ■ ^^^^ 
The differential {dV/dE^) and total decay width are singular at the lower limit of the 
photon energy, and this singularity which is present only in the | Mf^"' | contribution is 
due to the soft photon emission from charged lepton line. On the other hand, | M2 \ and 
Re [Mi^^'^M]] terms are free from this singularity. In this limit the Be TVrl decay 
can not be distinguished from the Be rPr decay. Therefore, in order to obtain a finite 
result for the decay width, we must consider both decays together. The infrared singularity 
arising in the | M"'^"' f contribution must be canceled with 0{a) virtual correction to the 
Be — > rUr decay. In the SM this cancellation explicitly was shown in [13]. 

In this work, the Be — > rPrl process is not considered as a 0{a) correction to the 
Be — > TPr decay, but rather a separate decay channel with hard photon radiation. There- 
fore we impose a cut off value on the photon energy, which will set an experimental limit on 
the minimum detectable photon energy. We consider the case for which the photon energy 
threshold is larger than 50, MeV i.e., E^ > amB^, where a > 0.01. An integration over all 
the possible values of the lepton energy E2 gives the total decay width as a function of the 
photon energy: 



r = 



647r2 



Be 



Vr 



cb 



c 



2HDM 



i: 



dx 



x{l — x) 



4 + 8r - 4r + lOx — lArx 
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+ Ar'^x - + 7rx^ + 3x^ + {1 - x){2 - 2r^ - 3x + rx + 2x'^)in{ 



1 — X, 



AfB. C^^^^M dx-^ (l-r-x) (Mx)x + Mx)(l + r - 2x)) 



- {l-x){h{x)x + h{x){2r-x)ynC—^) 

1 /■i-'^ , r, „ . . , „ . . ,91 1 
3)5 



+ ^ / dx 



x'^{2 + r - 2x)(l - r - x)^ >, 



(25) 



{l-xf ) 
where x = 2E^/mB^is the dimensionless photon energy, r = m1/m\^ and 5 = 2a. 

3 NUMERICAL ANALYSIS 

To calculate the decay width, explicit forms of the form factors /i and /2 are needed. These 
form factors were calculated in the framework of the light- front quark model in [14], and 
in the light-cone QCD sum rules [3], where in [3] it was found that, the best agreement is 
achieved by the following pole forms for the form factors: 

/i(0) , , 2^ /2(0) 



1 — /m\ ' 



1 — p^/ml ' 



where 



/i(0) = 0.44 ±0.04 GeV , 



43.1 GeV^ , 



(26) 



/2(0) = 0.21 ± 0.02 GeV , ml = 48.0 GeV^ . 

On the other hand, in evaluating the decay width, we have used the following set of param- 
eters: Gf = 1.17 .10-^ Gey-^ a = 1/137, rUb = A.%GeV, = lAGeV, tub, = &.2,GeV, 
rrir = 1.78 GeV, fs, = 0.36 GeV [15], Vcb = 0.04 [16] and, r(5c) = 0.46 x lO^^^ ^ ^^j^ 

In this regard we should also recall that the free parameters of the 2HDM model namely 
tanP, and Mh are not arbitrary, but there are some semiquantitative restrictions on them 
using the existing experimental data. The most direct bound on the charged Higgs boson 
mass comes from top quark decays, which yield the bound Mh > lATGeV for large tan(3 
[17]. For pure type-II 2HDM's one finds Mh > 300GeV^, coming from the virtual Higgs 
boson contributions to 6 — > 57 [18]. Furthermore, there are no experimental upper bounds 
on the mass of the charged Higgs boson, but one generally expects to have Mh < ITeV 



in order that perturbation theory remain vaUd [19]. For large tanP the most stringent 
constraints on tanP and Mh are actually on their ratio, tanP/Mn- The current limits come 
from the measured branching ratio for the inclusive decay B — > XtP, giving tanP/Mff < 
0.4:QGeV~^ [20], and from the upper limit on the branching ratio for B — > ri^, giving 
tanp/Mn < 0.38 GeV'^ [21]. 

For illustrative purposes we consider four values of tanP, namely tan(3 —5, 10, 30 and 
50 and let mH± — 150 GeV. Then we consider two values of mi/±, namely M^ff±=200, 400 
GeV, and we allow tanj3 to range between to 60. The results of this numerical analysis 
are graphically shown in figures 2-4. In these figures the differences between the 2HDM's 
and the SM are shown for two different fixed cut off values, i.e., S — 0.016 and 5 — 0.032 
both for Model I and Model II. 

The results for the differential decay branching ratio dBR(i?c — >■ Ti"^)/dx as a function 
ol X — 2E^/mB^ for different values of tanl3, Mh — 150 GeV are presented in Fig.2, while 
the branching ratio (BR) for B TP7 decay is shown in figure 3-4 as a function of Mh for 
various values of tan/3, and as a function of ton/3 for different values of Mh- Results are 
shown for Model I, and Mode II. 

It is observed that Model II gives both a bigger differential decay branching ratio, and 
a bigger branching ratio than the SM rates of (up to three orders of magnitude [Mh — 
150 GeV]) for large values of tanf3 > 20, while for small values of tan(3 < 10 results 
approaches its SM value. 

In model I the situation is somewhat totally different. Curves overlap with the SM results all 
the way. This behavior obviously reflects the if^ fermion couplings, which are proportional 

to cotjS in this model. 

In conclusion, this study shows that the branching ratios for B^. tuj could be at the 
level of 10~^ in the 2HDM, which may be detectable at the ongoing LHC. When enough Be 
events are collected, this decay will be able to provide alternative channel to extract new 
restrictions for the free parameters tanf3 and Mh+ of the 2HDM model. 
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Figure Captions 



Figure 1. : The relevant Feynman diagrams, responsible for Be tV^ decay. 
Figure 2. : The dependence of the differential Branching ratio of dBR(Sc — > Tv^)/dx as a 
function of the photon energy x = 2E^/mB^ for both models Model 1, and Model II. 
Figure 3. : The dependence of the Branching ratio on the charged Higgs boson mass at 
different values of tanl3 for both models Model I, and Model II. 

Figure 4. : The dependence of the Branching ratio on tan/3 at different values of the charged 
Higgs boson mass for both models Model I, and Model II. 
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